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1. Introduction 



Classical statistical physics deals with statistical systems in equilibrium. The ensemble theory 
offers a useful framework that allows to characterize and to work out the properties of this 
type of systems Q]]. Two fundamental distributions to describe situations in equilibrium are 
the Boltzmann-Gibbs (exponential) distribution and the Maxwellian (Gaussian) distribution. 
The first one represents the distribution of the energy states of a system and the second 
one fits the distribution of velocities in an ideal gas. They can be explained from different 
perspectives. In the physics of equilibrium, they are usually obtained from the principle of 
maximum entropy [2]. In the physics out of equilibrium, there have recently been proposed 
two nonlinear models that explain the decay of any initial distribution to these asymptotic 
equilibria dH. 

In this chapter, these distributions are alternatively obtained from a geometrical interpretation 
of different multi-agent systems evolving in phase space under the hypothesis of 
equiprobability. Concretely, an economic context is used to illustrate this derivation. Thus, 
from a macroscopic point of view, we consider that markets have an intrinsic stochastic 
ingredient as a consequence of the interaction of an undetermined ensemble of agents that 
trade and perform an undetermined number of commercial transactions at each moment. A 
kind of models considering this unknowledge associated to markets are the gas-like models 
0,0]. These random models interpret economic exchanges of money between agents similarly 
to collisions in a gas where particles share their energy. In order to explain the two before 
mentioned statistical behaviors, the Boltzmann-Gibbs and Maxwellian distributions, we will 
not suppose any type of interaction between the agents. The geometrical constraints and the 
hypothesis of equiprobability will be enough to explain these distributions in a situation of 
statistical equilibrium. 

Thus, the Boltzmann-Gibbs (exponential) distribution is derived in Section [2] from the 
geometrical properties of the volume of an N-dimensional pyramid or from the properties 
of the surface of an N-dimensional hyperplane @.@]- In both cases, the motivation will be a 
multi-agent economic system with an open or closed economy, respectively. The Maxwellian 
(Gaussian) distribution is derived in Section [3] from geometrical arguments over the volume 



2 



Stochastic Control 



or the surface of an N-sphere @,@]- In this case, the motivation will be a multi-particle gas 
system in contact with a heat reservoir (non-isolated or open system) or with a fixed energy 
(isolated or closed system), respectively. And finally, in Section |U the general equilibrium 
distribution for a set of many identical interacting agents obeying a global additive constraint 
is also derived |7]. This distribution will be related with the Gamma-like distributions found 
in several multi-agent economic models. Other two geometrical collateral results, namely the 
formula for the volume of high-dimensional symmetrical bodies and an alternative image of 
the canonical ensemble, are proposed in Section[5] And last Section embodies the conclusions. 



2. Derivation of the Boltzmann-Gibbs distribution 

We proceed to derive here the Boltzmann-Gibbs distribution in two different physical 
situations with an economic inspiration. The first one considers an ensemble of economic 
agents that share a variable amount of money (open systems) and the second one deals with 
the conservative case where the total wealth is fixed (closed systems). 



2.1 Multi-agent economic open systems 

Here we assume N agents, each one with coordinate x,-, i = 1, . . . , N, with x,- > representing 
the wealth or money of the agent i, and a total available amount of money E: 

Xi+x 2 -\ h x N _! + x N < E. (1) 

Under random or deterministic evolution rules for the exchanging of money among agents, 
let us suppose that this system evolves in the interior of the N-dimensional pyramid given 
by Eq. |(T). The role of a heat reservoir, that in this model supplies money instead of energy, 
could be played by the state or by the bank system in western societies. The formula for the 
volume Vjv(E) of an equilateral N-dimensional pyramid formed by N + 1 vertices linked by 
N perpendicular sides of length E is 

V N (E) = (2) 

We suppose that each point on the N-dimensional pyramid is equiprobable, then the 
probability f(xj)dxj of finding the agent i with money Xj is proportional to the volume formed 
by all the points into the (N — 1) -dimensional pyramid having the /th-coordinate equal to X;. 
We show now that /(*;) is the Boltzmann factor (or the Maxwell-Boltzmann distribution), 
with the normalization condition 

rE f{xi)dxi = l. (3) 



If the z'th agent has coordinate x,-, the N — 1 remaining agents share, at most, the money E — x,- 
on the (N — 1) -dimensional pyramid 

X\ + x 2 h + x i+1 h x N < E - x if (4) 

whose volume is V N _i (E — x,- ) . It can be easily proved that 

%(£) = /"%_!(£- xi)dxi. (5) 
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Hence, the volume of the N-dimensional pyramid for which the ith coordinate is between x, 
and Xj + dxj is V N _^(E — Xj)dx{. We normalize it to satisfy Eq. @, and obtain 

f{Xi) = V N (E) > (6) 



r.\N-l 
1 / - ' 



whose final form, after some calculation is 

f( Xi ) = NE~'[1 - | J , (7) 
If we call e the mean wealth per agent, E = Ne, then in the limit of large N we have 

li m (l-±i) ~e~ Xi/e . (8) 



N>i V E 

The Boltzmann factor e~ Xt ' e is found when N ^> 1 but, even for small N, it can be a good 
approximation for agents with low wealth. After substituting Eq. {§} into Eq. (7}, we obtain 
the Maxwell-Boltzmann distribution in the asymptotic regime N — > oo (which also implies 
E -S- oo): 

f(x)dx = X - e- x/e dx, (9) 

where the index i has been removed because the distribution is the same for each agent, and 
thus the wealth distribution can be obtained by averaging over all the agents. This distribution 
has been found to fit the real distribution of incomes in western societies 

This geometrical image of the volume-based statistical ensemble 1 7] allows us to recover the 
same result than that obtained from the microcanonical ensemble |8] that we show in the next 
section. 



2.2 Multi-agent economic closed systems 

Here, we derive the Boltzmann-Gibbs distribution by considering the system in isolation, that 
is, a closed economy. Without loss of generality, let us assume N interacting economic agents, 
each one with coordinate x,-, i = 1, . . . ,N, with x,- > 0, and where x, represents an amount of 
money. If we suppose now that the total amount of money E is conserved, 

x\ + x 2 H h x N _! + x N = E, (10) 

then this isolated system evolves on the positive part of an equilateral N-hyperplane. The 
surface area Sn(E) of an equilateral N-hyperplane of side E is given by 



Different rules, deterministic or random, for the exchange of money between agents can be 
given. Depending on these rules, the system can visit the N-hyperplane in an equiprobable 
manner or not. If the ergodic hypothesis is assumed, each point on the N-hyperplane is 
equiprobable. Then the probability /(x,)dx,' of finding agent i with money x, is proportional 
to the surface area formed by all the points on the N-hyperplane having the z'th-coordinate 
equal to Xj. We show that /(x,) is the Boltzmann-Gibbs distribution (the Boltzmann factor), 
with the normalization condition {5). 
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If the z'th agent has coordinate x- u the N — 1 remaining agents share the money E — Xj on the 
(N — l)-hyperplane 

x\ + x 2 h + x i+ i h x N = E - x if (12) 

whose surface area is Sjv-i(£ — *,-). If we define the coordinate 9^ as satisfying 



sinf?N= V^T' (13) 

it can be easily shown that 

fE Ay. 

S N (E) = S N ^(E-x,)—j-. (14) 
Jo sin fjv 

Hence, the surface area of the N-hyperplane for which the ith coordinate is between x,- and 
Xj + dxj is proportional to Sn_i(E — / sinS^j. If we take into account the normalization 

condition @, we obtain 

f(x) - 1 Su-rjE-Xj) 
f[X,) ~ S N (E) sin0 N ' 
whose form after some calculation is 

/(x,) = (N-l)£- 1 (l-^) , (16) 

If we call e the mean wealth per agent, £ = Ne, then in the limit of large N we have 

/ r-\N-2 

lim(l- = ) ~ e - x ' /£ . (17) 



N»i V E 

As in the former section, the Boltzmann factor e~ Xt ' £ is found when N»l but, even for small 
N, it can be a good approximation for agents with low wealth. After substituting Eq. (TJ into 
Eq. Jl6b , we obtain the Boltzmann distribution l(9) in the limit N — > oo (which also implies 
£ — > oo). This asymptotic result reproduces the distribution of real economic data fioll and 
also the results obtained in several models of economic agents with deterministic, random or 
chaotic exchange interactions I6l. ll2l.ll3l . 

Depending on the physical situation, the mean wealth per agent e takes different expressions 
and interpretations. For instance, we can calculate the dependence of e on the temperature, 
which in the microcanonical ensemble is defined by the derivative of the entropy with respect 
to the energy. The entropy can be written as S = — kN J °° f(x) ln/(x) dx, where f(x) is given 
by Eq. ^9) and k is Boltzmann's constant. If we recall that e = E/N, we obtain 

S(E) =fcNln-j| +kN. (18) 

The calculation of the temperature T gives 

T-i=( d 4) = k 4 = k -. (19) 



dE 



Thus e = kT, and the Boltzmann-Gibbs distribution is obtained in its usual form: 

f(x)dx = i e- x/kT dx. (20) 
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3. Derivation of the Maxwellian distribution 

We proceed to derive here the Maxwellian distribution in two different physical situations 
with inspiration in the theory of ideal gases. The first one considers an ideal gas with a variable 
energy (open systems) and the second one deals with the case of a gas with a fixed energy 
(closed systems). 



3.1 Multi-particle open systems 

Let us suppose a one-dimensional ideal gas of N non-identical classical particles with masses 
m{, with i = 1,...,N, and total maximum energy E. If particle i has a momentum m,o,-, we 
define a kinetic energy: 

K { = p? = im,p?, (21) 

where p,- is the square root of the kinetic energy K{. If the total maximum energy is defined as 
E = R 2 , we have 

V\ + Pi + ■ ■ ■ + Pn-i + P 2 n< R 2 - (22) 
We see that the system has accessible states with different energy, which can be supplied by 
a heat reservoir. These states are all those enclosed into the volume of the N-sphere given by 
Eq. {22). The formula for the volume Vjq(R) of an N-sphere of radius R is 
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V N (R) = — ^ R N , (23) 

where !(■) is the gamma function. If we suppose that each point into the N-sphere is 
equiprobable, then the probability /(p,)dp, of finding the particle i with coordinate p, (energy 
p 2 ) is proportional to the volume formed by all the points on the N-sphere having the 
ith-coordinate equal to p,-. We proceed to show that /(p,) is the Maxwellian distribution, 
with the normalization condition 

_ R f( Pi )d Pi = l. (24) 

If the z'th particle has coordinate p,-, the (N — 1) remaining particles share an energy less than 
the maximum energy R 2 — p? on the (N — 1) -sphere 

Pi + v\ ■ ■ ■ + Pi-i + Pm • • • + p 2 n < R 2 - pI (25) 



whose volume is V^j- 1 ( y R 2 ~ p})- It can be easily proved that 

V N (R) = j\ N _ l (^R 2 -p 2 )dp l . (26) 
Hence, the volume of the N-sphere for which the ith coordinate is between p,- and p,- + dpi is 



Vn— 1 ( \/R 2 ~ P 2 )dpi- We normalize it to satisfy Eq. l|24l , and obtain 



ft 



- v N{R) > W 
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whose final form, after some calculation is 



v 2 \ — 

/(p,) = C N R- 1 (l-^) 2 , (28) 



with 

7 nY{ 

For N>1, Stirling's approximation can be applied to Eq. ( f29t , leading to 



Cn = — (29) 



1 N 

limC-^T. (30) 
If we call e the mean energy per particle, E = R 2 = Ne, then in the limit of large N we have 

N-l 

lim I 1 "!! ) 2 ^ p ' /2e - (3D 

n»i I R 2 J 

The factor e~P 2 / 2i? is found when N>1 but, even for small N, it can be a good approximation 
for particles with low energies. After substituting Eqs. ll3"0"t-{3Tl into Eq. {28}, we obtain the 
Maxwellian distribution in the asymptotic regime N — > co (which also implies E — > co): 



f(p)dp=^^- e e-P- /2€ dp, (32) 

where the index i has been removed because the distribution is the same for each particle, and 
thus the velocity distribution can be obtained by averaging over all the particles. 

This newly shows that the geometrical image of the volume-based statistical ensemble Cj 
allows us to recover the same result than that obtained from the microcanonical ensemble Oj] 
that it is presented in the next section. 



3.2 Multi-particle closed systems 

We start by assuming a one-dimensional ideal gas of N non-identical classical particles with 
masses nij, with i = 1, . . . , N, and total energy E. If particle i has a momentum m,w,, newly 
we define a kinetic energy Kj given by Eq. {21), where p,- is the square root of Kj. If the total 
energy is defined as E = R 2 , we have 

V\ + Vl + " ' + Vn-i + Pn = R 2 - (33) 

We see that the isolated system evolves on the surface of an N-sphere. The formula for the 
surface area Sjv(-R) of an N-sphere of radius R is 

N 

Sn(R) = ^w~R N -\ (34) 

where T(-) is the gamma function. If the ergodic hypothesis is assumed, that is, each point 
on the N-sphere is equiprobable, then the probability f(fi)dpj of finding the particle i with 
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coordinate p; (energy p?) is proportional to the surface area formed by all the points on the 
N-sphere having the ith-coordinate equal to p,-. Our objective is to show that /(p;) is the 
Maxwellian distribution, with the normalization condition i24l . 

If the z'th particle has coordinate p,-, the (N — 1) remaining particles share the energy R 2 — p? 
on the (N — 1) -sphere 

Pi + V\ ■ ■ ■ + Pi-i + rf+i ■ ■ ■ + P 2 n = R 2 ~ Pl (35) 



whose surface area is Sj[_i(^/R 2 — p?). If we define the coordinate 6 as satisfying 

R 2 cos 2 6 = R 2 — p 2 , (36) 



then 



It can be easily proved that 



Rd8 = q . (37) 

Pi v ' 
R 2 ■ 



(l-SV/2 



rn/2 

S N (R) = / S N ^ 1 (Rcos9)Rde. (38) 

./-7T/2 

Hence, the surface area of the N-sphere for which the ith coordinate is between p,- and p,- + dpi 
is Sj V f_i(J?cos6)Rd0. We rewrite the surface area as a function of p z -, normalize it to satisfy 
Eq. l|24t , and obtain 

! S N _ 1 (jR 2 -p 2 ) 
/Pi = V , — - , (39) 

SnW (1-^)1/2 

whose final form, after some calculation is 

„2 N-3 



with 



/(p,)=C N R- 1 (l-|L) 2 , (40) 

_ i r(f) 

CN -7^rTV)- (41) 

For N>1, Stirling's approximation can be applied to Eq. d4Tb , leading to 

i fW 

lim C N - ^\ T - (42) 
If we call e the mean energy per particle, E = R 2 = Ne, then in the limit of large N we have 

lim | 2 ~e-rf /2e . (43) 

n»i ^ R 2 y v ; 

As in the former section, the Boltzmann factor e~P7/ 2e is found when N S> 1 but, even for 
small N, it can be a good approximation for particles with low energies. After substituting 
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Eqs. d2)-(43} into Eq. (40), we obtain the Maxwellian distribution (32) in the asymptotic 
regime N — > oo (which also implies E — > oo). 

Depending on the physical situation the mean energy per particle e takes different 
expressions. For an isolated one-dimensional gas we can calculate the dependence of e on the 
temperature, which in the microcanonical ensemble is defined by differentiating the entropy 
with respect to the energy. The entropy can be written as S = —kNj^ /(p)ln/(p) dp, where 
/(p) is given by Eq. (32) and k is the Boltzmann constant. If we recall that e = E/N, we obtain 

S(E) = hcNln ( 1) + ifcN(ln(2/r) - 1). (44) 



2 V N . 
The calculation of the temperature T gives 

Thus e = kJ '/2, consistent with the equipartition theorem. If p 2 is replaced by jtnv 2 , the 
Maxwellian distribution is a function of particle velocity, as it is usually given in the literature: 



(v)dv = J^e- m ^^ T dv. (46) 



4. General derivation of the equilibrium distribution 

In this section, we are interested in the same problem above presented but in a general way. 
We address this question in the volume-based statistical framework. 

Let b be a positive real constant (cases b = 1, 2 have been indicated in the former sections). If 
we have a set of positive variables (xi,X2, ■ ■ ■ , xjv) verifying the constraint 

x\ + x\ + ■ ■ ■ + + x b N < E (47) 

with an adequate mechanism assuring the equiprobability of all the possible states 
(x\,Xi, . . . , xjv) into the volume given by expression (47), will we have for the generic variable 
x the distribution 

f(x)dx~e- 1/b e- xb/be dx, (48) 

when we average over the ensemble in the limit N, E — > oo, with E = Ne, and constant e?. 
Now it is shown that the answer is affirmative. Similarly, we claim that if the weak inequality 
(47) is transformed in equality the result will be the same, as it has been proved for the cases 
b = l,2inRefs. EH]. 

From the cases b = 1,2, (see Eqs. (6) and (27)), we can extrapolate the general formula 
that will give us the statistical behavior f(x) of the generic variable x, when the system runs 
equiprobably into the volume defined by a constraint of type l l47t . The probability f(x)dx of 
finding an agent with generic coordinate x is proportional to the volume Vjv_i((E — x 6 ) 1 / 6 ) 
formed by all the points into the (N — 1) -dimensional symmetrical body limited by the 
constraint (E — x b ). Thus, the N-dimensional volume can be written as 

V N (E 1/b )= [ EUl V N _ 1 ((E-x b ) 1/fc )dx. (49) 
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Taking into account the normalization condition J f(x)dx = 1, the expression for f(x) is 
obtained: 

The N-dimensional volume, V^(b,p), of a fo-symmetrical body with side of length p is 
proportional to the term p N and to a coefficient gb(N) that depends on N: 

V N (b,p)=g b (N)p N . (51) 

The parameter b indicates the original equation {47} that defines the boundaries of the volume 
Vjv (b, p). Thus, for instance, from Eq. {2), we have g b —i(N) = 1/Nl. 

Coming back to Eq. {50}, we can manipulate V^((E — x 11 ) 1 ^) to obtain (the index b is omitted 
in the formule of V^): 

/ \ - 

V N {{E-x b ) 1 ' b )= gb {N) [(E-x^Y = gb (N)ET fl-yj • 02) 

If we suppose E = Ne, then e represents the mean value of x b in the collectivity, that is, 
e =< x b >. If N tends toward infinity, it results: 

N 

x ■ -zVbe 



lim 1 - — = e- x /oe . (53) 
N»i \ E ' 

Thus, 

V N ((E - x b ) 1/b ) = V N (E 1/b ) e- x " /be . (54) 

Substituting this last expression in formula {50}, the exact form for f(x) is found in the 
thermodynamic limit (N, E — > oo): 

f(x)dx = c b £- 1/b e- x " /b£ dx, (55) 

with c b given by 

Cb = toW-V. (56) 

g b (N)N^ 

Hence, the conjecture {48} is proved. 

Doing a thermodinamical simile, we can calculate the dependence of e on the temperature 
by differentiating the entropy with respect to the energy. The entropy can be written as S = 
— kN Jq°° f(x) ln/(x) dx, where f(x) is given by Eq. {55} and k is the Boltzmann constant. If we 
recall that e = E/N, we obtain 

S(E) = ^ln(|)+^(l-fcInc fc ), (57) 

where it has been used that e =< x b > = x b f(x)dx. 
The calculation of the temperature T gives 
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Thus e = kT/b, a result that recovers the theorem of equipartition of energy for the quadratic 
case b = 2. The distribution for all b is finally obtained: 

f(x)dx = c b (JlV /b e-*"' a dx. (59) 



4.1 General relationship between geometry and economic gas models 

If we perform the change of variables y = e^^x in the normalization condition of f(x), 
Io° fi x )d x = 1/ where f(x) is expressed in {55}, we find that 



(60) 



If we introduce the new variable z = y 'lb, the distribution fix) as function of z reads: 



/(z)dz = zs" 1 e~ 2 dz. (61) 
Let us observe that the Gamma function appears in the normalization condition, 

b 1 



[ f( z )dz=-^r [ z5- 1 e - 2 dz = ^ T r(M =1. (62) 

./o JW ./o m-I U; 

This implies that 



Cfc = —/TV ( 63 ) 

r(i) 

By using Mathematica the positive constant is plotted versus b in Fig. 1. We see that 
lim^o c b — °°/ an d that limj^K, C{, = 1. The minimum of C;, is reached for b = 3.1605, taking 
the value c b = 0.7762. Still further, we can calculate from Eq. (63[ the asymptotic dependence 
of Cj, onb: 

Hm Cfc =/X v W/»(l-A + ...), (64) 
lim c b = b- l/b (l + l + ...), (65) 



b 

where 7 is the Euler constant, 7 = 0.5772. The asymptotic function i64t is obtained 
after substituting in (63) the value of r(l/fo) by (1/b — 1)!, and performing the Stirling 
approximation on this last expression, knowing that 1/b — > 00. The function (651 is found 
after looking for the first Taylor expansion terms of the Gamma function around the origin 
x = 0. They can be derived from the Euler's reflection formula, r(x)T(l — x) = 71/ sin(7rx). 
We obtain Y(x — ► 0) = + T'(l) + ■ ■ ■ . From here, recalling that r'(l) = —7, we get 
r(l/b) = b — 7+ • • • , when b — > 00. Although this last term of the Taylor expansion, —7, 
is negligible we maintain it in expression J65). The only minimum of c b is reached for the 
solution b = 3.1605 of the equation ip(l/b) + logb + b — 1 = 0, where xp(-) is the digamma 
function (see Fig. 1). 
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Figure 1. Normalization constant cj, versus b, calculated from Eq. l63l . The asymptotic behavior is: 
limt_jo c b = 00 , and lim^oo cj, = 1. This last asymptote is represented by the dotted line. The minimum 
of Cj is reached for b = 3.1605, taking the value = 0.7762. 




Let us now recall two interesting statistical economic models that display a statistical behavior 
given by distributions nearly to the form ioTl , that is, the standard Gamma distributions with 
shape parameter 1/b, 

f(z)dz = -L-zT-'e- z dz. (66) 

ECONOMIC MODEL A: The first one is the saving propensity model introduced by 
Chakraborti and Chakrabarti [11]. In this model a set of N economic agents, having each 
agent i (with i = 1, 2, ■ ■ ■ , N) an amount of money, «,-, exchanges it under random binary (2, /) 
interactions, («,, Uj) — > (u'j, «'■), by the following the exchange rule: 

u\ = Xui + e(l - A)(u,- + uj), (67) 
u'j = Aiij + e(l - A) (uj + Uj), (68) 

with e = (1 — e), and e a random number in the interval (0, 1). The parameter A, with < A < 
1, is fixed, and represents the fraction of money saved before carrying out the transaction. Let 
us observe that money is conserved, i.e., m, + Uj = !(' + u'y hence in this model the economy 
is closed. Defining the parameter n(A) as 

n W = (69) 



and scaling the wealth of the agents as z = nu/ < u >, with < u > representing the average 
money over the ensemble of agents, it is found that the asymptotic wealth distribution in this 



12 



Stochastic Control 



system is nearly obeying the standard Gamma distribution IT3.[T^1 

f(z)dz = z"- 1 e- 2 dz. (70) 

The case n = 1, which means a null saving propensity, A = 0, recovers the model of 
Dragulescu and Yakovenko 1 10] in which the Gibbs distribution is observed. If we compare 
Eqs. (701 and {66), a close relationship between this economic model and the geometrical 
problem solved in the former section can be established. It is enough to make 

n = 1/b, (71) 
z = z, (72) 

to have two equivalent systems. This means that, from Eq. (7Tt , we can calculate b from the 
saving parameter A with the formula 

(73 > 

As A takes its values in the interval (0, 1), then the parameter b also runs in the same interval 
(0, 1). On the other hand, recalling that z = x b /be, we can get the equivalent variable x from 
Eq. (72), 

e iVi 



[ 



< u > 

b : 



(74) 



where e is a free parameter that determines the mean value of x in the equivalent geometrical 
system. Formula ((74) means to perform the change of variables w, — > x,-, with i = 1,2, • • ■ , N, 
for all the particles /agents of the ensemble. Then, we conjecture that the economic system 
represented by the generic pair (A, u), when it is transformed in the geometrical system given 
by the generic pair (b, x), as indicated by the rules \73\ and <74t , runs in an equiprobable form 
on the surface defined by the relationship (47), where the inequality has been transformed in 
equality. This last detail is due to the fact the economic system is closed, and then it conserves 
the total money, whose equivalent quantity in the geometrical problem is E. If the economic 
system were open, with an upper limit in the wealth, then the transformed system would 
evolve in an equiprobable way over the volume defined by the inequality <47t . although its 
statistical behavior would continue to be the same as it has been proved for the cases b = 1, 2 
inRefs. fill]. 

ECONOMIC MODEL B: The second one is a model introduced in Q3]. In this model a set 
of N economic agents, having each agent i (with i = 1, 2, ■ • ■ , N) an amount of money, w,, 
exchanges it under random binary (i,f) interactions, («,-, Uj) — > («', up, by the following the 
exchange rule: 

u ■ = Uj - An, (75) 
u'j = Uj + Au, (76) 

where 

Am = r](x.j — Xj) ecvxi — [1 — jy(x; — xy)] ewxj, (77) 

with e a continuous uniform random number in the interval (0,1). When this variable is 
transformed in a Bernouilli variable, i.e. a discrete uniform random variable taking on the 
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values or 1, we have the model studied by Angle 11711 . that gives very different asymptotic 
results. The exchange parameter, to, represents the maximum fraction of wealth lost by one 
of the two interacting agents (0 < CO < 1). Whether the agent who is going to loose part of 
the money is the z'-th or the j-th agent, depends nonlinearly on (x, — Xj), and this is decided 
by the random dichotomous function 17(f): t](t > 0) = 1 (with additional probability 1/2) 
and f](t < 0) = (with additional probability 1/2). Hence, when x,- > Xj, the value i] = 1 
produces a wealth transfer from agent i to agent with probability 1/2, and when x,- < Xj, the 
value f] = produces a wealth transfer from agent to agent i with probability 1/2. Defining 
in this case the parameter n(co) as 

and scaling the wealth of the agents as z = nu/ < u >, with < 11 > representing the average 
money over the ensemble of agents, it is found that the asympt otic wealth distribution in this 
system is nearly to fit the standard Gamma distribution 



isympto 



f(z)dz = j^-y z"- 1 e" 2 dz. (79) 

The case n = 1, which means an exchange parameter co = 3/4, recovers the model of 
Dragulescu and Yakovenko fioll in which the Gibbs distribution is observed. If we compare 
Eqs. {79} and {66}, a close relationship between this economic model and the geometrical 
problem solved in the last section can be established. It is enough to make 

n = l/b, (80) 
z = z, (81) 

to have two equivalent systems. This means that, from Eq. J80t . we can calculate b from the 
exchange parameter to with the formula 

As co takes its values in the interval (0, 1), then the parameter b runs in the interval (0, 2). It 
is curious to observe that in this model the interval co 6 (3/4, 1) maps on b 6 (1,2), a fact 
that does not occur in MODEL A. On the other hand, recalling that z = x b /be, we can get the 
equivalent variable x from Eq. (81}, 

l/b 

(83) 



L < u > J 



where e is a free parameter that determines the mean value of x° in the equivalent geometrical 
system. Formula ((83} means to perform the change of variables m, — > x,-, with i = 1,2, • • ■ , N, 
for all the particles /agents of the ensemble. Then, also in this case, we conjecture that 
the economic system represented by the generic pair (A, u), when it is transformed in the 
geometrical system given by the generic pair (b, x), as indicated by the rules (82) and \83\ , runs 
in an equiprobable form on the surface defined by the relationship (47), where the inequality 
has been transformed in equality. As explained above, this last detail is due to the fact the 
economic system is closed, and then it conserves the total money, whose equivalent quantity 
in the geometrical problem is E. If the economic system were open, with an upper limit in the 
wealth, then the transformed system would evolve in an equiprobable way over the volume 
defined by the inequality (47), although its statistical behavior would continue to be the same 
as it has been proved for the cases b — 1, 2 in Refs. 0||^]. 
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5. Other additional geometrical questions 

As two collateral results, we address two additional problems in this section. The first one 
presents the finding of the general formula for the volume of a high-dimensional symmetrical 
body and the second one offers an alternative presentation of the canonical ensemble. 



5.1 Formula for the volume of a high-dimensional body 

We are concerned now with the asymptotic formula (N — > oo) for the volume of the 
N-dimensional symmetrical body enclosed by the surface 

x\ + x\ + ■ ■ ■ + x b N _ r + x b N = E. (84) 

The linear dimension p of this volume, i.e., the length of one of its sides verifies p ~ E 1 ^. As 
argued in Eq. d5lt , the N-dimensional volume, V^/(b,p), is proportional to the term p N and to 
a coefficient gt,(N) that depends on N. Thus, 

V N (b,p)= 8b (N)p N , (85) 

where the characteristic b indicates the particular boundary given by equation {84). 
For instance, from Equation ((2), we can write in a formal way: 

ft-xCN) - f ^ Ty (86) 

From Eq. {23}, if we take the diameter, p = 2R, as the linear dimension of the N-sphere, we 
obtain: 

N 

gb=2(N) = v (87) 

These expressions {86} and {87} suggest a possible general formula for the factor gb(N), let us 
say 

gb (N) = 



f + 1 



where a is a fe-dependent constant to be determined. For example, a = 1 for b = 1 and a = n/i 
for b = 2. 



In order to find the dependence of a on the parameter b, the regime N — > oo is supposed. 
Applying Stirling approximation for the factorial ( y ) ! in the denominator of expression {88} , 
and inserting it in expression {56}, it is straightforward to find out the relationship: 

c b = (ab)- 1/b . (89) 



From here and formula {63}, we get: 

'1 



b +1 



(90) 
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Figure 2. The factor g b (N) versus b for N = 10,40, 100, calculated from Eq. j9lj. Observe that 
g b (N) = for b = 0, and \im b ^ g b (N) = 1. 




that recovers the exact results for b = 1,2. The behavior of a is monotonous decreasing when 
b is varied from b = 0, where a diverges as a ~ l/b + ■ ■ ■ , up to the limit 6 — > oo, where a 
decays asymptotically toward the value «oo = = 0.5614. 

Hence, the formula for gb{N) is obtained: 

r(f + i) 

It would be also possible to multiply this last expression d91~b by a general polynomial K(N) in 
the variable N, and all the derivation done from Eq. d88t would continue to be correct. We omit 
this possibility in our calculations. For a fixed N, we have that gi, (N) increases monotonously 
from gt,{N) = 0, for b = 0, up to gb{N) = 1, in the limit b — > oo (see Fig. 2). For a fixed b, we 
have that gb(N) decreases monotonously from gt(N) = 1, for N = 1, up to gb(N) = 0, in the 
limit N->oo (see Fig. 3). 

The final result, that has been shown to be valid for any for the volume of an 

N-dimensional symmetrical body of characteristic b given by the boundary ( f84t reads: 



r(i + i 



v N (b, P ) = y; r \ P N , (92) 



with p ~ £ 1/fc . 
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Figure 3. The factor g b (N) versus N for b = 10,40, 100, calculated from Eq. j9lj. Observe that 
g b (N) = 1 for N = 1, and lim^^, g h (N) = 0. 




5.2 A microcanonical image of the canonical ensemble 

From Section[2] here we offer a different image of the usual presentation that can be found in 
the literature U about the canonical ensemble. 

Let us suppose that a system with mean energy E, and in thermal equilibrium with a heat 
reservoir, is observed during a very long period t of time. Let E; be the energy of the system 
at time i. Then we have: 

Ej + E 2 + ■ ■ ■ + E T _\ + E T = t ■ E. (93) 

If we repeat this process of observation a huge number (toward infinity) of times, the different 
vectors of measurements, (Ej, E2, . . . , E T _\, E T ), with < E, < t ■ E, will finish by covering 
equiprobably the whole surface of the r-dimensional hyperplane given by Eq. {93}. If it is 
now taken the limit t — > 00, the asymptotic probability p(E) of finding the system with an 
energy E (where the index i has been removed), 

p(E) ~ e- E/£ , (94) 

is found by means of the geometrical arguments exposed in Section [2] Doing a 

thermodynamic simile, the temperature T can also be calculated. It is obtained that 

E = kT. (95) 

The stamp of the canonical ensemble, namely, the Boltzmann factor, 

p{E) ~ e- E/kT , (96) 

is finally recovered from this new image of the canonical ensemble. 
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6. Conclusion 

In summary, this work has presented a straightforward geometrical argument that in a certain 
way recalls us the equivalence between the canonical and the microcanonical ensembles in 
the thermodynamic limit for the particular context of physical sciences. In the more general 
context of homogeneous multi-agent systems, we conclude by highlighting the statistical 
equivalence of the volume-based and surface-based calculations in this type of systems. 

Thus, we have shown that the Boltzmann factor or the Maxwellian distribution describe the 
general statistical behavior of each small part of a multi-component system in equilibrium 
whose components or parts are given by a set of random linear or quadratic variables, 
respectively, that satisfy an additive constraint, in the form of a conservation law (closed 
systems) or in the form of an upper limit (open systems), and that reach the equiprobability 
when they decay to equilibrium. 

Let us remark that these calculations do not need the knowledge of the exact or microscopic 
randomization mechanisms of the multi-agent system in order to attain the equiprobability. 
In some cases, it can be reached by random forces jfj, in other cases by chaotic Il3l Il9ll or 
deterministic causes. Evidently, the proof that these mechanisms generate equiprobability 
is not a trivial task and it remains as a typical challenge in this kind of problems. 

The derivation of the equilibrium distribution for open systems in a general context has also 
been presented by considering a general multi-agent system verifying an additive constraint. 
Its statistical behavior has been derived from geometrical arguments. Thus, the Maxwellian 
and the Boltzmann-Gibbs distributions are particular cases of this type of systems. Also, 
other multi-agent economy models, such as the Dragalescu and Yakovenko's model the 
Chakraborti and Chakrabarti's model JTHl and the modified Angle's model IT3l . show similar 
statistical behaviors than our general geometrical system. This fact has fostered our particular 
geometrical interpretation of all those models. 

We hope that this framework can be useful to establish other possible relationships between 
the statistics of multi-agent systems and the geometry associated to such systems in 
equilibrium. 
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